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The structures of the configuration space of the six-vertex models with 
various boundaries and boundary conditions are investigated, and it is de- 
rived that the free energies depend on the boundary conditions, and that 
they are classified by the fractal structures. The " n-equivalences" of the 
boundary conditions are defined with a property that the models with n- 
equivalent boundary conditions result in the identical free energy. The con- 
figurations which satisfy the six-vertex restriction are classified, through the 
n-equivalences, into sets of configurations called islands. It is derived that 
each island shows a fractal structure when it is mapped to the real axis. Each 
free energy is expressed by the weighted fractal dimension (multi-fractal di- 
mension) of the island. The fractal dimension of the island has a strict 
relation with the maximum eigenvalue of the corresponding block element of 
the transfer matrix. It is also found that the fractal dimensions of islands 
take the absolute maximum when the boundary is "alternative"; in this case, 
the free energy equals to the solution obtained by Lieb and Sutherland. 
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1 Introduction 



Recently, Korepin and Zinn- Justin obtainedfl] the free energy of a six-vertex 
model with a specific boundary condition. Their free energy is written in 
terms of elementary functions and different from that obtained, by Lieb[2- 
4] and Sutherland [5], with the condition that the boundary is cyclic. This 
seems unusual compared to other systems, such as the Ising models with 
finite range interactions, where the free energy per site is independent of the 
boundary condition in the thermodynamic limit. 

This fact can be understood by the following argument. The six- vertex 
configurations have a one-to-one correspondence with configurations of con- 
tinuous lines on the lattice. The boundary configuration determines the 
number of lines on the lattice, and configurations with different line numbers 
are disjoint [6] each other in the configuration space. Let "island" be the set of 
all the configurations we can find with a specific condition on the boundary. 
The partition function is obtained as the sum of Boltzmann weights on the 
island. Now it seems natural that models with different boundary conditions 
and hence defined on different islands result in different free energies. 

In this paper, we investigate the story shortly explained above, especially 
the relation between the free energy and the structures of the set of config- 
urations. It is derived that the free energies are classified according to the 
islands, each island has a fractal structure and its fractal dimension has a 
strict relation to its free energy. 

First we introduce new equivalence class of boundary conditions called the 
" n-equivalence" . Two boundary conditions yield just the same free energy 
if they are n-equivalent for some finite n. (This argument is general and 
independent of Hamiltonians.) 

Next a specific boundary condition so-called the alternative boundary 
condition is introduced and it is proved that this boundary condition yields 
the free energy /ls which was obtained, by Lieb and Sutherland, with the 
cyclic boundary condition. In Theorem 1, it is derived that the free energy 
is identical with / LS in any six- vertex models on arbitraly domain D, where 
the boundary is an closed and continuous curve, provided that the boundary 
condition is alternative. This result also shows that the free energy is still 
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intensive even if the boundary effect remains relevant in the thermodynamic 
limit. 

Next we investigate the structures of the configuration space of the six- 
vertex models. It is derived that each island shows generally non-integral 
Hausdorff dimension when it is appropriately mapped to the real axes, i.e. 
they are fractal. Each island results in its own fractal dimension. There exists 
an simple relation between the free energy and the weighted fractal dimension 
(multifractal dimension) defined on this fractal island. This structure yields 
the fact that the free energy is dependent on the boundary condition. We 
can also show that the fractal dimensions of islands take the maximum when 
the boundary condition is alternative. The free energy from this island is 
/ls- The free energies of six-vertex models with more complex boundary 
conditions can also be classified from the structure of the configuration space. 

It should be noted that the configuration space is simply the set of con- 
figurations, between which we can find relations coming from the six-vertex 
restriction. The fractal dimension appears when the island is mapped to the 
real axis referring the relations between the configurations. We simply use 
the term the fractal dimension of the island. Each fractal dimension d of an 
island corresponds to the high temperature limit of the maximum eigenvalue 
A of an irreducible block element of the corresponding transfer matrix, and 
the larger d corresponds to the larger A. 

Let us consider the square lattice and assign an arrow on each bond. The 
arrows are arranged so that two arrows come in and the other two go out 
at each site. Then there exist six types of possible local arrow arrangements 
as shown in Fig.l. In this paper we are going to use the term 'vertex' as a 
site and four bonds around it. Each vertex is assumed to have finite energy. 
The energies are assumed to be unchanged by reversing all the arrows on 
the lattice. Then we have three energy parameters and hence three types of 
Boltzmann weights a, b and c assigned to each vertex (see again Fig.l). We 
also introduce j3 = 1/ksT, k-Q is the Boltzmann constant, T is the tempera- 
ture and N is the total number of sites. The partition function is 

N 

E IK*s (!) 

config. i=l 

where q is the energy of the i-th vertex and ^config. is taken over all the 
possible configurations. The free energy is —f3f = linijv^co iV -1 log Z. 
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One can assign a line on the bonds with arrows going down or going 
left (Fig.l). The six- vertex restricton corresponds to the restricton that the 
lines don't intersect each other and continue till they reach to the boundary 
of the lattice. Namely, each vertex configuration just corresponds to a line 
configuration of continuous lines on the lattice. 

We consider a domain D that is finite and simply-connected. The bound- 
ary of D is a closed and continuous line j(t) = (x(t),y(t)) (0 < t < 1) which 
satisfies 7(^1) 7^ l(t 2 ) if h 7^ t 2 except 7(0) = 7(1). We also assume that the 
sites are on the points (niai,n 2 a 2 ), where n\ and n 2 are integers, a\ and a 2 
are the lattice spacings. The vertices belong to D when the corresponding 
sites lie inside of D. The bonds/vertices in D are called as the boundary 
bonds/vertices, when they have non-zero intersection with 7. The sites of 
the boundary vertices are called as the boundary sites. These definitions are 
rather tedious but necessary because we will take the following limit: we fix 
the line 7 and take the thermodynamic limit a l7 a 2 — > 0. This corresponds 
to taking the limit w, h — > 00 where w and h are the number of columns and 
rows in D. We also assume that the number of vertices on the boundary 
and the energy contributions from the boundary, respectively, divided by the 
total number of sites vanish in the thermodynamic limit. 



§ 2 The n-equivalence 

First we introduce the n-equivalence. Let us consider the set of sites which lie 
in D and can be reached from the boundary sites by just n steps (n bonds) at 
minimum; we call these sites as the n-boundary sites. Furthermore consider 
the set of bonds between (n — 1)- and n-boundary sites, and call them as 
the n-boundary bonds. Let {Tj} be the set of all the possible configurations 
on the n-boundary with a fixed configuration T on the boundary of D. Two 
boundary conditions Y and V are called n-equivalent when {Tj} = {T^} as a 
set of n-boundary configurations. 

The n-equivalence is a generalization of the concept of boundary con- 
dition. The definition is general and can be introduced to models other 
than the six-vertex model. For example, all the boundary configurations 
are 1-equivalent to each other in the Ising model because its spin states are 
independent of their nearest-neighbors. 

Proposition 1 Suppose that the boundary conditions Y and Y' are n-equivalent 
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with some finite n. Then the two free energies with F and V are identical in 
the thermodynamic limit. 

Proof: Let us write the partition function as Z = Y^iBiZi. The factor Z t 
is the partition function of the set of vertices inside the n-boundary with 
fixed ry The factor Bi is the contribution from outside the n-boundary, 
with fixed n-boundary configuration iy The factor BiZi is the partition 
function of the system with fixed ry From the assumption (written in the 
last paragraph of introduction) that the boundary energy is irrelevant in 
the thermodynamic limit, together with the fact that n is finite, we obtain 
log Zi = —j3Nfi + o(N) and log-Bj = o(N). The index i runs from 1 to i max 
where i max is the number of possible configurations on the n-boundary with 
fixed boundary configuration Y. We have i max < 0(2^ ) where N' is the 
number of vertices on the n-boundary. Then we obtain 

^o g Z = llogQTi^) 

i 

1 mi R 7 * max R 

= ^o g B l z l[ i + z 2ff+ E §-e-w<-n+«*] 

i=2 1 1 i=mi+l 1 

- h (N^oo), (2) 

where f t = f x (1 < i < mi), /i < /< (mi + 1 < i). ■ 



§ 3 The free energies on the domain D with 
fixed boundary configurations 

In this section, we derive that the free energy of the six-vertex model on the 
domain D is /ls if the boundary line configuration of D is a specific one 
called "alternative". The main result of this section is Theorem 1. 

First let us consider a rectangle R with w columns and h rows. Assume 
that h and w are even. Let us consider a row of vertices in R. Let r\ = 
{x±, . . . , x n } and 77' = {x'i, . . . , x' n } be line configurations with n lines (0 < 
n < w) on the upper and the lower edges of the row, respectively. The 
symbol {x±, X2, ■ ■ ■ , x n } denotes that we find lines (find down arrows) just on 
the Xk-th bond (k = 1, . . . ,n), and we don't find lines (find up arrows) on 
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the x-th bond with x ^ %k- The (rf, ^-element of the transfer matrix V is 
defined as 




w 



(3) 



config. fc=l 



where (k — 1, . . . , w) is the energy of the fc-th vertex on the row and here 
the sum I] con fi g . is taken over all the possible line configurations on the row 
with fixed 77 and rj'. 

If we identify the bond at the right end with that at the left end in each 
row, which we call the cyclic boundary condition in the horizontal direction, 
the transfer matrix is the same for all rows. If the boundary condition is also 
cyclic in the vertical direction the partition function Z is written as 



where Aj are the eigenvalues of V and Ai > |Aj| for all i. 

In reference [7], for example, one can find a nice review of the six- vertex 
model (ice-type model) on R, with the cyclic boundary condition in both 
horizontal and vertical directions, together with the review on the transfer 
matrix method applied to it. Following [7] set A = (a 2 + b 2 — c 2 )/2ab. The 
transfer matrix is block diagonalized according to the number of lines which 
intersect each row. The maximum eigenvalue of the transfer matrix with the 
cyclic boundary lies in the block element with or w lines (A > 1), with 
w/2 lines (A < 1). The free energy / for the case A > 1 is / = min(e a ,eb), 
where exp(— j3e a ) = a and exp(— /3e&) = b. Hence we concentrate on the case 
A < 1. From the symmetry of the system, we can conclude that the number 
of lines which intersect a column should be also h/2. 

Lemma 1 Suppose that the boundary condition of R is cyclic in the hori- 
zontal direction. Then all the line configurations with n lines on the upper 
boundary row (and also on the lower boundary row) of R are 2n-equivalent 
to each other. 

Proof: Let {xi,X2, ■ ■ ■ ,x n }, X{ < Xi+i, be a line configuration on the upper 
boundary row of R. When we have a line arrangement {xi, X2, x$, x±, x$} = 



Z = trV h = J2 X i ~ X i ( h 



00), 



(4) 
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{2,3,5,6,9}, just as shown in Fig. 2(a), one can introduce the shift of lines 
{2, 3, 5, 6, 9} -> {1, 2, 3, 5, 6} -> {1, 2, 3, 4, 5}. When we generally have 
{1, 2, . . . , i, Xj, Xj+i, . . . , x n _fc}, i + 1 < a;j, after the /c-th step, one can in- 
troduce the line arrangement {1, 2, . . . , i, i + 1, x^, . . . , x n _( fe+1 )} as the next 
one. We need n — I steps for the shift {x±, x 2 , ■ ■ ■ , £„} — > {1, 2, . . . , n} if 
xi < n < xi+i, and this means {xi, ^2, . . . , x„} is at most n-equivalent to 
{1,2, ... ,n}. Then arbitrary line configurations {x[, x 2 , . . . ,x' n } are again 
possible on the 2n-boundary. | 

Let us consider a specific boundary configuration where the arrows going 
right and left, or going up and down, appear alternatively on the boundary. 
We assume that the arrow on the first row (column) goes left (down). This ar- 
row configuration corresponds to the line configuration in which lines appear 
on every other boundary bonds. When we fix the boundary configuration as 
alternative we call it as the " alternative boundary condition" . 

Lemma 2 Suppose that the boundary condition of R is fixed as alternative 
on the right and the left edges. Then all the line configurations with n lines 
on the upper (and also on the lower) boundary of R are An- equivalent to each 
other 

Proof: We need 2n steps for the shift {x±, x 2 , ■ ■ ■ , x n } — 2, . . . , n} when we 
introduce the shift {1,2,..., i, Xj,Xj + i, . . . , x n -k}^ {1,2, ... ,i, i+l, Xj, . . . , x n ^(k+i)} 
and the "parallel shift" {x[, x' 2 , . . . , x' n }— >{x' n , x[, . . . , x' n _ x }, as shown in Fig. 
2(b), every other horizontal rows. We also need 2n steps at most for the shift 
to an arbitrary configuration {1,2,..., n}—>{x", x 2 ', . . . , x^} when we intro- 

(k + 1)} and the "identical shift" {x[, x' 2 , . . . , x' n }— >{x[, x' 2 , . . . , x' n } every 
other horizontal rows. Then arbitrary line configurations {x'[, x 2 , . . . , x'^} 
are again possible on the 4n-boundary. | 

We assume the existence of the free energy in the limit h — > oo w — > oo 
and also in the limit w — > oo h — > oo for any fixed boundary conditions. 
Then we obtain the following result. 

Proposition 2 Suppose that the boundary condition of R is cyclic in the 
horizontal direction, or otherwise suppose that the boundary line configura- 
tions are fixed as alternative on the right and the left edges of R. Then the 
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free energy of the six-vertex model on the lattice R is unique in the limit 
(w,h) — > (00,00), i.e. the limit is independent of the order of two limiting 
procedures w — > 00 and h — > 00, and also independent of the ratio w/h. 

Proof: Suppose first that the boundary condition is cyclic in the horizontal 
direction. The elements of the transfer matrix V are non-negative because 
they are sums of Boltzmann weights. Let V n be the block element of V for a 
fixed line number n. All the line configurations with n lines are 2n-equivalent 
to each other, so all the elements of V^ n are positive. That means V% n is 
irreducible, and hence V n is irreducible. Then the Frobenius theorem works 
and we know the following. Let \i(w) be the i-th eigenvalue of V n and let 
Ai(tf) > \{w). Then we have \\{w) > and Xi(w) > \\i(w)\. We also 
know that all the elements of the eigenvector for the eigenvalue Xi(w) are 
positive, i.e. the projections satisfy (xi, . . . ,x n |max)> where \x±, . . . ,x n ) 
is the state with the fixed upper line arrangement {x±, . . . , x n } and |max) is 
an eigenstate for the maximum eigenvalue of V n . 

The partition function Z with finite w and h is written as 

Z = (2^1, • • • , x n \ V \x± : ... : x n ) 

= c 1 X 1 (w) h + J2cMw) h (5) 

where {xi, . . . , x n } and {x[, . . . , x' n } are the line configurations on the upper 

and the lower edges of R, respectively, and c± — (x[, . . . , x' n \ma.x) (max|a;i, . . . , x n ) > 

0. We will consider the limit of 

■^-logZ = J-bgci + — logAi(tw) + \z'(h,w) 
hw liw w h 

z'(h,w) = log[l + E-(^)1 1/w (6) 

i>2 C l A l 

The factor z'(h, w) is a decreasing function of h. Taking the limit h — > 00 
first and then w — > 00 we find the convergence of the factor w~ 1 logAi(w). 
Taking the limit w — > 00 first we find from (6) that the factor z'(h,w) is 
convergent, and hence anyway bounded by some constant. Then the con- 
vergence of the free energy in (6) is uniform as h — > 00, hence the limit is 
independent of the order of the limiting procedures w — > 00 and h — > 00, and 
also independent of the ratio w/h when we take the limit with < w/h < 00. 
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If the boundary condition is alternative on the right and the left edges 
of R, we can introduce a transfer matrix V = V1V2 where V\ and V2 are 
the transfer matrix of successive two rows of R, respectively. The partition 
function on R is written as a linear combination of the terms \i(w) h ^ 2 , where 
Xi(w) are now the eigenvalues of V — ViV 2 , and we obtain the same result. | 

Note that the n-equivalence corresponds to the irreducibility of the block 
element V n : line configurations are n-equivalent to each other if they belong 
to the same bases set of an irreducible block element of V. This is a key to 
generalize the results for the cyclic boundary into our general cases. 

It is obvious that the proof can be generalized to the cases with arbitrary 
periodic line configurations on the right and the left edges and boundary 
conditions which is n-equivalent to it with fixed n. 

When the boundary is cyclic in the horizontal direction \±(w) is known 
and we have lim,^^ w" 1 log Xi(w) = — /3/ls- Hence the next proposition is 
actually proved in the proof of proposition 2. However here we again express 
the proof using only the n-equivalence and the uniform convergence of the 
free energy. 

Proposition 3 Suppose that the boundary condition of R is cyclic in the 
horizontal direction, and the line configurations on the upper and the lower 
edges are fixed with w/2 lines. Then the free energy of the six-vertex model 
on the lattice R is /ls, in the limit (w, h) — > (00, 00). 

Proof: The line configuration on the upper edge is w-equivalent to the ar- 
bitrary line configurations in the block element with w/2 lines, and hence 
w-equivalent to the cyclic boundary with w/2 lines. We can take the limit 
h — > 00 with fixed w, and we obtain the free energy with (w, h) — > (w, 00). 
The resulted functions are the same for the two alternative and the cyclic 
boundary conditions, on the upper and the lower edges. Next taking the 
limit w — > 00, we finally find the free energy which is equal to f LS in both 
cases. 1 

In the proof of the next proposition we don't need the explicit form of 
AiH- 

Proposition 4 Suppose that the four boundaries of R are fixed as alterna- 
tive. Then the free energy of the six-vertex model on the lattice R is /ls, in 
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the limit (w, h) — > (oo, oo). 

Proof: Let us identify the left edge, which is fixed as alternative, with the 
right edge, which is also fixed as alternative. Then we obtain a tube of a 
height h with the upper and the lower edges as its boundaries. 

The alternative boundaries on the upper and the lower edges are 4w- 
equivalent to the cyclic boundary with w/2 lines. Taking the limit h — > oo 
with fixed w, we obtain the free energy which is identical to that obtained 
with the cyclic boundary in the vertical direction and with fixed alternative 
boundaries on the left and right edges. Next taking the limit w — > oo, we 
again obtain / LS . | 

Now we can obtain a rather general sufficient condition to have /ls- Let 
D be the domain defined in introduction. We introduce, on the boundary, 
the line density p which is the number of lines over the number of the bonds. 

Theorem 1 The free energy of the six-vertex model on the lattice on D is 
fhs, if the line density p on the boundary equals to 1/2. 

Proof: We first assume that the boundary of D is convex. Let D be a 
rectangle of a width w and a height mh satisfying D C D$. The rectangle 
.Do can be divided into m small rectangles Ri (i = 1, . . . , m) of the height h; 
the lower edge of Ri coincide with the upper edge of Ri + \. Let D' = U™ =1 i?^, 
where the R\ are the rectangles of the height h and of the width w i: satisfying 
R'i C Ri and R\ C D, and satisfying that the number of columns of each R' i: 
Wi = riiv, takes its maximum with the restriction that R[ C D. The sum of 
the line numbers on the upper and the right-hand sides of R[ equals to the 
sum of the line numbers on the lower and the left-hand sides, because of the 
line conservation property. 

Fix the line density on the boundary of D' being p = 1/2. Then there 
exist line configurations where the line density on the edges of each R\ are 
1/2. We are able to find the maximum eigenvalue of the transfer matrices of 
each R[ when the boundary line density of R[ is 1/2. 

Only the maximum eigenvalue of the transfer matrices remains when we 
take the limit 02 — > (which corresponds to h — > 00) with fixed 7, and next 
taking the limit ai — > (w — > 00), each rectangle R^ yields the free energy 
/ls as shown in Proposition 4. Then we obtain fjy = /ls where fo 1 is the 
free energy on the lattice D' . 
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The result is independent of the ratio h/w and the ratio can be taken 
sufficiently small. (Note that the line density of both D and D' are 1/2.) 
The ratio of the energy contribution, the contribution from D\D' over that 
from D, goes to zero because the line 7, which determines the boundary of 
D, is continuous. Hence we obtain \Jd~ fwl < e for arbitrary positive e. The 
region D can be divided into sums of small regions with convex boundaries. 1 

This means that the free energy is still additive, i.e. proportional to the 
volume of the system, in the situation where the boundary condition remains 
relevant inside D even in the thermodynamic limit. 

It should be noted that we didn't need diagonalize successive transfer 
matrices but it was sufficient for us to consider the dominant configulations 
which are n-equivalent to each other. 

Here we fixed the line density p to be 1/2 but our argument is valid for 
any < p < 1. 

§ 4 Fractal structure hidden in the configu- 
ration space 

Next we consider the structure of the configuration space W and its relation 
to the free energy of the system. 

Let us consider the set of vertices on the domain D. Each bond of the 
vertices can take two states, right or left (or up or down). Let k be the 
number of bonds, then 2 k configurations are generally possible. We call this 
set of configurations as the configuration space W . We finally take the limit 
k — > 00. 

Let u be a subset of bonds in D and Y be an arrow configuration on u. 
Let U U (Y) be the set of all the possible configurations in W with fixed u and 
T. Obviously U U '(T') C U U (T) when u C u' and f = T on u. 

Next let us introduce the set Wi which is the set of all the possible con- 
figurations with a fixed arrow configuration or a specific condition on the 
boundary of D. We call Wi as an island. 

The space W is simply the set of configurations. We will introduce a map 
F : W 1— * R and introduce F(Wi) = which is a set of real numbers coming 
from the island Wi. We show that Ei generally have non- integral Hausdorff 
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dimension in the thermodynamic limit. 

We are going to find a mapping which has the following properties: 

• Each configuration has equal contribution to the Hausdorff measure. 

• Adding a new vertex, Ei get the next detailed structure and then the 
asymptotic behavior of the free energy in the thermodynamic limit 
k — > oo determines the fractal dimension of E { . 

In this paper we introduce the following rule as a possible mapping F to 
the real axis. Fix an order of bonds and let the j-th bond be <jj and assign 
Sj to each <jj where Sj = +1 (or 0) if the arrow on <x,- goes down/left (or 
up/right). Let x = Y$=i s j2~- ? and we find a mapping F from configurations 
in W to real numbers x G [0, 1]. 

Our choice works especially in the rectangle with a fixed width, as shown 
in Lemma 3, and has the following advantages: 

• Each Ei for the rectangle is a generalization of the Cantor set, which 
is one of the most typical fractal sets. 

• Some special but interesting cases of E^ become a family of graph- 
directed sets, and hence one can use the general results obtained for 
them. 

The second point can considerably simplifies arguments. We will, for exam- 
ple, estimate the Hausdorff dimension g?h and the box-counting dimension 
g?b- To estimate g?h is usually extremely difficult but we have, in the graph- 
directed sets, g?h — d-Q and d-Q can be obtained much easier than d^. 

We are now considering the measure of E iy especially for the cases where 
D is a rectangle. Counting the measure of E^ corresponds to calculating the 
partition function with a specific boundary condition in the high temperature 
limit, because we have Z = Sconfig. II e - ^ — > Sconfig. 1 (P — V^b^ - ► 0). 
The partition function at finite temperatures is obtained through assigning 
a finite weight to each configuration and integrating over the islands. 

Let us consider the first k bonds u = {a±, <t 2 , • • • , &k\ and fix the arrow 
configuration on it as Y = (si, s 2 , . . . , Sfc). Then we find that U U (T) is an 
interval in [0, 1] with the volume 

|C4(r)| = sup{\x-x'\\x,x' E U u (T)} = 2- k . (7) 
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Let If ' = [{I - l)/2 fc , l/2 k ] (I = 1, . . . , 2 k ). Each If > corresponds to a fixed 
arrow configuration on the first k bonds. We concentrate on the intervals 
1^ which satisfy x G 1^ for some x G E i} and write them as iff , I G 
{li, l 2) . . . ,l mk } C {1, 2, 3, ... , 2 k }. Here m k is the number of such intervals. 
Let us introduce = \JY=i 4r • From now on we omit the index i for the 

(k) (k) 

intervals and simply write I- / = // for convenience. 

Lemma 3 Let R(w) be a rectangle with a fixed width w. Then the corre- 
sponding and 1^ satisfy 
i; if) D £(*+i> 

ii) The locations of small intervals I^ 2 ^ in each 1^ are determined by s k 
and s k -2w, where k = w + (h! - l)(2w + 1) + (2w' - 1), h',w' G N, 1 < w' < w. 

Proof: Let us introduce an order of bonds as shown in Fig.3, the vertical w 
bonds on the first row as ci, (x 2 , . . . , a w from right to left, the other 2w + 1 
bonds of the vertices of the first row as cr w +i, o- w+2 , ■ ■ ■ , o~3w+i again from 
right to left, and vise versa. 

Let us consider, for example, the vertex with the bonds o"n, a 2 i, o 22 
and (723. Fix (sn,S2i) = (0, 1), then two line arrangements are possible on 
a 22 and a 23 as (522,^23) = (0,1) or (1,0), because only the six arrow ar- 
rangements are permitted as shown in Fig.l. Hence ij 21 ^ fl E^ consists of 
two intervals of the width 1/2 23 beginning from a\ + 1/2 23 and a\ + 2/2 23 , 
where a\ is the left end of corresponding I\ . Only the following finite 
number of patterns are possible, corresponding to the six types of vertices, 
such as (s k - 2w ,s k ,s k+1 ,s k+2 ) =(0,0,0,0), (0,1,0,1), (0,1,1,0), (1,0,0,1), 
(1, 0, 1, 0) and (1, 1, 1, 1), where k = w + (h'- l)(2w + 1) + (2w' - 1), h! G N, 
w' — 1, 2, 3, . . . , w. Then we obtain i) and ii) provided that the width w is a 
constant. 1 

Together with the fact that linn^oo = 0, the property i) certifies the 
convergence of E^ (k — > 00) on the real axis. 

(k) 

Each of these ^ corresponds to a line configuration satisfying the six- 
vertex restriction on the first k bonds. The number of configurations m k is 
obtained from the high temperature limit of the partition function with fixed 
k. We know the asymptotic behavior of m k (for example from the argument 
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in [7]) as 



m k = const, x rf + o(rf) (N — > oo), 
n = lim Xi{w) 1/w (8) 

T — >oo 

where Aj(w) is the maximum eigenvalue corresponding to the line density pi, 
N is the number of sites, 2N = k except the irrelevant contributions from 
the boundary. 

Proposition 5 Let Wi be an island for R(w), pi be the line density on the 
boundary of R{w) and Ei = F(Wi). The box-counting dimension of Ei is 
d B = logn/2 log 2. 

Proof: The box-counting dimension is defined through the volume of the in- 
tervals as rf = (l/2 k )~ dB and we have d B = logr;/21og2. | 

When we consider the six-vertex model on arbitraly domain D, Lemma 
3 i) is not necessarily valid throughout the thermodynamic limit and Ei 
generally has more complicated structures on the real axis. However we 
know from Theorem 1 that the free energy does not depend on the shape of 
D and anyway obtained from the corresponding \\{w). Here we would like 
to note that n -> (4/3) 3/2 at p = p x = 1/2 and d B -> log(4/3) 3/2 /2 log 2 
when w — > oo. 

We can also estimate the Hausdorff dimension of E{. The s-dimensional 
Hausdorff measure of E t is defined as 

H s (E t ) = hmmfJ2\U 3 \ s , (9) 

^ j 

where {Uj} is a collection of subsets satisfying E { C U^Uj, sup\Uj \ = 5 > 0, 
and the infimum is over all {Uj}. The Hausdorff measure 7i s {E,j) is an 
decreasing function of s. It jumps from oo to at a specific value of s = g?h — 
dimH-Ei which we call the Hausdorff dimension of Ei. In some nice cases the 
Hausdorff measure takes finite and non-zero values at s = du- The Hausdorff 
dimension is countably stable: diniH-E = dimn Ei = sup 1< j <00 dimH-Ei, 
which means we generally know that the islands with the maximum Hausdorff 
dimension determines dim^E even in the thermodynamic limit where we 
generally have infinite number of islands. 
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Proposition 6 The Hausdorff dimension of Ei is du = logrj/21og2. 

Proof: Let us consider an 1^ and write 1^ fl E\ k+V> = {///}. If we have a 
real number s which satisfies 



m 



'\I' V \ S = \rf k) \ s , m is the number of //, in l[ k \ (10) 



for all k and /, we obtain H s (Ei) — 1, — HkLo^i^ ( see f° r example [8], 
Th.1.15), which means the Hausdorff dimension of E,- L is s. Especially if the 
intervals 1^ are uniformly spaced, the condition (10) is written in terms of 
the number of intervals m& as 

m kl +k 2 , 1 y = ,J_y , X 

that is 

m fcl+fc2 = 2 fc2S m fcl = 2^ +k >>. (12) 

The proof in the Th.1.15 can be straightforwardly generalized to our case 
where (I) the intervals //, satisfy (8) instead of (12) and (II) the intervals 
are not equally spaced but number of small intervals in a larger interval J 
is almost proportional (statistically proportional) to its width |J|. Then the 
remaining point is to show the property (II): 

Let R(w, h) be the set of vertices on the first h rows of R(w). Let k\ and 
k 2 be the number of bonds in R(w, hi) and R(w, h 2 ), respectively. Assume 
hi < h 2 and hence ki < k 2 . We will take the limit h — > oo with fixed hi and 
h 2 . 

We know that the number of intervals l\ kl ^ and l\ k2 ^ behaves ~ 
r 2kl + o(r 2kl ) and rrik 2 — r 2k ' 2 + o(r 2fc2 ), respectively, when ki and k 2 are 
sufficiently large. The domain R(w , h 2 )\R(w , hi) is again an rectangle with 
the hight h 2 — h\. Then the number of intervals in an interval lj kl \ 
= l/2 fc i, behaves r 2 ^~ k ^ + (r 2 ^' k ^). Let re be the number of l\ k2) 
in J = \j Q /2 k \ (jo + Aj)/2 fel ] = bo/2 fel , (jo + l)/2 fel ] U • • • U [(j + Aj - 
l)/2 kl , (jo + Aj)/2 kl ]. The relation between \J\ — Aj/2 kl and re is almost 
linear when k 2 is sufficiently large. | 

We assume the alternative boundary condition on the right and the left 
edges of R(w). In this case the result is obtained from a more general argu- 
ment. The system becomes a family of graph-directed sets defined by a set 
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of IFS (Iterated Function System). The Hausdorff dimension g?h of a family 
of graph-directed sets generally coincides with their box-counting dimension 
o?b (see for example [10]). Hence we don't need estimate du but it is sufficient 
for us to calculate d B . 

Proposition 7 With the condition above, {Ei} with fixed number of lines 
form a family of graph- directed sets. 

Proof: The rectangle R(w) is composed by iterations of a unit which is a 
small rectangle with w columns and two rows. The set of possible configura- 
tions on the lower edge of a unit is determined by the configuration on the 
upper edge of the same unit. Corresponding to this fact, we find that each 
Ei is obtained by iterations of set of contractions (IFS), where each IFS is 
determined by the configuration on the upper edge of each unit ( = the lower 
edge of the previous unit). The system also satisfies the transitivity condi- 
tion, which means arbitrary line configuration on an edge with n lines can 
be obtained from the other arbitrary line configuration on an previous edge 
with n lines after finite times of iterations of the unit, because all the line 
configurations with n lines are 4n-equivalent to each other. Taking h — > oo, 
we find that the Ej is a member of a family of graph-directed sets {Ej} 
(jn < j < jn+i) where each index j = j n ,j n + 1, . . . ,j n +i — 1 corresponds to 
a line configuration with n lines on the first edge of R(w). | 

The well known property of graph-directed sets, that the dimension of Ej 
are equal for j satisfying j n < j < j n+ ±, corresponds to the 4n-equivalences 
of the line configurations on the first row. This proposition can be straight- 
forwaldly generalized to the cases with other boundary conditions. Note that 
systems on a rectangle R(w) can be treated in terms of the graph-directed 
sets and the corresponding transfer matrix is of a finite degree. 

At last we consider the finite temperature case. Let us consider a general- 
ized fractal dimension so-called the multifractal dimension, which is defined 
as 

^ ) = fe /?-llog(l/2*) ' (13) 

where iuf is the Boltzmann weight coming from l\ k \ The free energy of the 
system is generally written as 

Pf(j3)/2]og2 = (J3-l)iM(J3). (14) 
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The sum in (13) is written as the sum over all the islands as 

E< = EE°' ) K°Y, (is) 
i j i 

where is the weight from the Z-th interval of the j-th island and J2i ^ 
is the sum with fixed j. The multifractal dimension Pj(/3) is introduced 
restricting the sum in (13) to the sum in the j-th island. We have J2i ^ 1 — 
0((l/2 k )~ dB 'i) (k — > oo) where dsj is the box-counting dimension of the j-th 
island, d B ,i > ^bj- We also have that the high temperature limit T — > oo 
(/3 — > 0) gives /i(0) = c?b,i = <fe- The dimension is obtained from 

the maximum eigenvalue of the corresponding transfer matrix. This is the 
relation between the structure of the configuration space and the free energy 
of the system. 

§ 5 Discussions and summary 

The six-vertex model is usually solved on a rectangle R in terms of the 
transfer matrix which is block diagonalized according to the number of lines. 
The free energy with a fixed number of lines is obtained from the maximum 
eigenvalue of a block element. 

Theorem 1, and its straightforward generalization to the cases with line 
density p (0 < p < 1), says that the six- vertex model on the domain D with 
fixed p on the boundary results in the free energy which is obtained from the 
maximum eigenvalue of the corresponding block element of the corresponding 
transfer matrix. When p = 1/2, we obtain /ls- 

The high temperature limit T — > oo ((3 — > 0) of the free energy gives 
the fractal dimension of the corresponding island. When we can find many 
islands corresponding to each boundary configuration, the island with the 
maximum fractal dimension determines the free energy of the whole system 
when A < 1 in sufficient high temperatures. 

If A > 1, we have two types of frozen phases [7], where two types of line 
configurations become dominant. In this case the free energy is constant and 
the generalization from a rectangle R to the domain D and all the other 
corresponding arguments given in this paper become trivial. 

The relation between the free energy and the structure of the configura- 
tion space is summarized in terms of the multifractal dimension. Each island 
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has its own multifractal dimension. In the limit T — > oo the multifractal 
dimension converges to the box-counting dimension. At low temperatures 
we possibly find transitions of phases. 

We can introduce boundary conditions, such as the cyclic boundary con- 
dition, in which various boundary configurations T are admitted. If our 
vertex energies and the temperature satisfy A < 1, and if we can find V 
which is n-equivalent to the alternative boundary for some finite n, then the 
free energy is /ls- If A > 1 and if T with and with w lines are admitted 
the system result in the trivial frozen phase. 

This result explains several results already obtained and also determines 
the exact free energies of six-vertex models which have not yet been solved, 
as follows. We have several exact calculations about six-vertex models with 
different conditions on the boundaries and they result in just the same free 
energy / LS . For example, Batchelor et al. solved[ll] (by the Bethe ansatz 
method) the six- vertex model with the periodic and "free" boundary condi- 
tions, respectively, in two directions. Owczarek and Baxter solved[12] (by the 
Yang-Baxter relation) the six-vertex model on a rectangular domain R with 
the periodic and the anti-periodic boundary conditions in the horizontal and 
the vertical directions, respectively. The anti-periodic boundary condition 
means that the arrow at one end of a row should be rightward (leftward) if 
that on another end of the same row is leftward (rightward). These results 
were obtained through a direct calculations and considerable efforts were 
needed. However we can find that the alternative boundary is realized with 
the restrictions both in [11] and [12], and hence this is now obvious from our 
result that the free energies of these systems are /ls- Furthermore we can 
easily construct a bunch of n-equivalent cases which have not yet been solved 
and are extremely difficult to solve directly by the Bethe ansatz method or 
the Yang-Baxter relation, but we know that the free energies of all of these 
cases should be / L s- 

We also have another example[l] where the alternative boundary is pro- 
hibited. This case yields a distinct free energy as just noted in introduction. 

Our argument is an extreme generalization of the fact that the maximum 
eigenvalue of the transfer matrix determines the free energy. We showed the 
correspondence between fractal dimensions of sets on the real axis and the 
maximum eigenvalues of block elements of the transfer matrix. Actually we 
have derived that the islands corresponding to the rectangle with a finite 
width become convergent figures on the real axis and show non-integral frac- 
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tal dimensions. We already know from Theorem 1 that the free energy is 
independent of the shape of the boundary of D. The free energy for D is ob- 
tained from the maximum eigenvalue of the corresponding block element of 
the corresponding transfer matrix. It is to be hoped that the results obtained 
for the rectangles will be generalized in some form to the cases of aribitrary 
domain D. 

At last it should be noted that we found a correspondence between a 
class of fractal figures and the transfer matrices and hence a class of lattice 
models. 

The basic structure studied in this paper is general. It is interesting to 
move to consider the self-similarity at the critical point or the equivalences 
of Hamiltonians. It is also interesting to consider the relations to dynamical 
systems. The properties of statistical systems should have more strict rela- 
tions to structures hidden in the configuration space. 

The author would like to thank Professor V. Korepin for discussions. 
The author would also like to thank Professor T. Hattori for stimulating 
comments and discussions. 
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Figure Captions 

Fig. 1 Six vertices, corresponding line configurations and their Boltz- 
mann weights. 

Fig. 2 Shift of lines introduced in the proof of Lemma 1 and Lemma 2. 
Fig. 3 The order of bonds in R(w), where the width is fixed as w = 5. 
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